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Abstract — This paper considers the problem of energy- 
efficient transmission in multi-flow multihop cooperative 
wireless networks. Although the performance gains of 
cooperative approaches are well known, the combinatorial 
nature of these schemes makes it difficult to design efficient 
polynomial-time algorithms for joint routing, scheduling 
and power control. This becomes more so when there is 
more than one flow in the network. It has been conjectured 
by many authors, in the literature, that the multiflow 
problem in cooperative networks is an NP-hard problem. 
In this paper, we formulate the problem, as a combinatorial 
optimization problem, for a general setting of fc-flows, and 
formally prove that the problem is not only NP-hard but 
it is o(n^/'^ inapproxmiable. To our knowledge*, these 
results provide the first such inapproxmiablity proof in 
the context of multiflow cooperative wireless networks. We 
further prove that for a special case of k = 1 the solution is 
a simple path, and devise a polynomial time algorithm for 
jointly optimizing routing, scheduling and power control. 
We then use this algorithm to estabUsh analytical upper 
and lower bounds for the optimal performance for the 
general case of k flows. Furthermore, we propose a 
polynomial time heuristic for calculating the solution for 
the general case and evaluate the performance of this 
heuristic under different channel conditions and against 
the analytical upper and lower bounds. 

I. Introduction 

In a wireless network, a transmit signal intended for 
one node is received not only by that node but also by 
other nodes. In a traditional point-to-point system, where 
there is only one intended recipient, this innate property 
of the wireless propagation channel can be a drawback, 
as the signal constitutes undesired interference in all 
nodes but the intended recipient. However, this effect 
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also implies that a packet can be transmitted to multiple 
nodes simultaneously without additional energy expen- 
diture. Exploiting this broadcast advantage, broadcast, 
multicast and multihop unicast systems can be designed 
to work cooperatively and thereby achieve potential 
performance gains. As such, cooperative transmission 
in wireless networks has attracted a lot of interest not 
only from the research community in recent years [1], 
[2], [4], [5], [6], [7], [8] but also from industry in the 
form of first practical cooperative mobile ad-hoc network 
systems [10]. The majority of the work in the cooper- 
ative literature has so far focused on the single flow 
problem, though recently there has been an increased 
interest in considering multiflow settings in cooperative 
networks [13], [14], [15], [16], [17], [9]. 

We consider a time-slotted system^ in which the nodes 
that have received and decoded the packet are allowed to 
re-transmit it in future slots. During reception, nodes add 
up the signal power (EA) received from multiple sources. 
Details of EA, and possible implementations have been 
extensively discussed in prior work [2], [4], [5], [14]. 
A key problem in such cooperative networks is routing 
and resource allocation, i.e., the question which nodes 
should participate in the transmission of data, and when, 
and with how much power, they should be transmitting. 
The problem is further complicated when there is more 
than one flow going through the network at the same 
time. 

We focus on the problem of minimum-energy multi- 
flow cooperative transmission in this paper, where there 
are k source-destination pairs, with each source node 
wanting to send a packet to its respective destination 
nodes, in a multihop wireless network. Other nodes in 
the network, that are neither a source nor a destina- 
tion, may act as relays to help pass on the message 
through multiple hops, provided they have already de- 
coded the message themselves and they are not trans- 
mitting/receiving any other messages at the time. The 
transmission is completed when all the destination nodes 
have successfully received their corresponding messages. 



Without loss of generality, we assume unit time-slots. Thus the 
terms energy and power are used interchangeably throughout this 
work. 
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It has been noted in the literature ([8], [18]) that a 
key tradeoff in cooperative settings is between the total 
energy consumption and the total delay measured in 
terms of the number of slots needed for all destination 
nodes in the network to receive the message. Therefore, 
we take delay into consideration and focus on the case 
where there is a delay constraint, whereby the destination 
node(s) should receive the message within some pre- 
specified delay constraint. We therefore formulate the 
problem of performing this transmission in such a way 
that the total transmission energy over all transmitting 
nodes is minimized, while meeting a desired delay 
constraint on the maximum number of slots that may be 
used to complete the transmission. The design variables 
in this problem determine which nodes should transmit, 
when, and with what power. 

We furthermore assume that the nodes are memory- 
less, i.e., accumulation at the receiver is restricted to 
transmissions from multiple nodes in the present time 
slot, while signals from previous time slots are discarded. 
This assumption is justified ([8], [18]) by the limited 
storage capability of nodes in ad-hoc networks, as well as 
the additional energy consumption nodes have to expand 
in order to stay in an active reception mode when they 
overhear weak signals in preceding time-slots. 

The main contribution of the work presented in this 
chapter is as follows: It has been conjectured in the 
literature that the problem of jointly computing schedules, 
routing, and power allocation for multiple flows in 
cooperative networks is NP-hard [16], [17], [15]. In this 
chapter we formulate the joint problem of scheduling, 
routing and power allocation in a multiflow cooperative 
network setting and formally prove that not only it is NP- 
hard, but it is also o{v}/'^~^) inapproximable. (i.e., unless 
P = NP, it is not possible to develop a polynomial time 
algorithm for this problem that can obtain a solution that 
is strictly better than a logarithmic-factor of the optimum 
in all cases). We are not aware of prior work on multiflow 
cooperative networks that shows such inapproximability 
results. We further prove that for a special case of A; = 1, 
the solution is a simple path and devise an optimal 
polynomial time algorithm for joint routing, scheduling 
and power control. We establish analytical upper and 
lower bounds based on this algorithm and propose a 
polynomial-time heuristic, the performance of which is 
evaluated against those bounds. 

The rest of this paper is organized as follows: In 
section II we provide a mathematical formulation of the 
problem. In section III we consider the special case of 
A; = 1 and prove the solution is a simple path and 
can be found optimally in polynomial time. The inap- 
proximabUty results are presented in section IV using 



reduction from minimum graph coloring problem. We 

establish analytical upper and lower bounds for optimal 
performance in section V. A polynomial-time heuristic is 
proposed in section VI and its performance is evaluated 
under different channel conditions and against the per- 
formance bounds. Concluding remarks are summarized 
in section VIII. 

II. Problem Formulation 

Consider a network, G, with a total of n nodes, 
/ = {1, ... n}. Assume we have r source nodes, labeled 
S = {si,S2, Sr}, and r corresponding destination 
nodes, V = {di,d2, ■■■,dr}. The source-destination 
nodes can be thought of as pairs, {(sfc, dk)}^.^^, all with 
the same delay constraint T. The goal is to deliver a 
unicast message from each source to its corresponding 
destination, possibly using other nodes in the network 
as relays. The objective is to do so using the minimum 
amount of sum transmit power and within the delay 
constraint. 

We consider a cooperative wireless setting with EA 
and consider signal-to-intereference-plus-noise (SINR) 
threshold model, [14], [2], [12], [11]. That is, in order 
for node i to be able to decode message k at time t, the 
following inequality needs to be satisfied: 

Here s/- (t) is the set of nodes transmitting the message k 
at time t, hij is a constant between and 1 representing 
the channel gain between node i and j, and N and 9 
are constants representing the noise and the decoding 
threshold respectively. 
Equation (1) can be re- written as 

n r n 

hjiP% ''-^Plt -0N>0, (2) 

j=l q=l u=l 

where pf^ is the power used by node i at time t to 
transmit message k. 

The system is memoryless, meaning although we are 
allowed to accumulate the same message from multiple 
sources during each time slot, we cannot accumulate over 
time. The relays are half-duplex, meaning they cannot 
transmit and receive simultaneously. The relays cannot 
transmit more than one message at the same time either. 

In order to apply ideas driven by the rich literature 
on multicommodity flows [19] to our problem, we need 
to somehow introduce the notion of delay constraint 
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into the multicommodity setting. What follows is a 

transformation of our network graph that would allow for 
the multicommodity flow technique to be applied, while 
observing the delay constraint: For a delay constraint 
T, map the given network to a layered graph with T 
layers as shown in Figure 1. Place a copy of all the 
nodes in the network on each of the layers. Connect 
each node, on each layer, to its corresponding copy on its 
neighboring layers with an edge weight of 0. Also create 
directed edges between each node, on each layer k, and 
the nodes on the next layer k + 1, with edge weights 
representing the amount of power required to transmit 
the message from the node on the top level to the node 
on the bottom level, as a whole. Notice that there is no 
edge between the nodes on the same level. Call the new 
graph G'. Assign the nodes corresponding to the source 
nodes of G on level 1 of G' as source nodes in G' and 
the destination nodes on level T of G', corresponding to 
destination nodes in G, as destinations in G', as shown in 
the figure. Similar transformations have been used in the 
Uterature in the context of multiflow transmission [16]. 



Without loss of generality, we assume unit length time 
slots. The nodes who want to transmit are to do so at 
the begiiming of each time slot, and the decoding (by 
nodes who receive enough information during that time 
slot) will happen by the end of that time slot. Let z^^ 
be an indicator binary variable that indicates whether or 
not node i decodes the message k during time slot t, as 
per inequality in equation (1). In other words, we define 
zff to be 1, if node i decodes message k during time 
slot t, and otherwise. Let p| be the transmit power 
used by node i at each time t to transmit message k. We 
define another binary variable that is 1 if node i is 
allowed to transmit message k at time t, and otherwise. 
A node is allowed to transmit during a particular time 
slot, if it has already decoded that message in previous 
time slots, and it's not receiving or transmitting any other 
messages during that time slot. Notice that being allowed 
to transmit does not necessarily mean that a transmission 
actually occurs. To take care of actual transmissions, let 
us define v^^ to be a binary variable that is 1 if node i 
transmits message k at time t, and otherwise. 



min Ptotal = Er=i Er=i ELi Pit 



(3) 



s.t. 



1. pf^ > 0, yi,t,k 

4. (-M)(l-4)<yf„ Vi,t 

5. p% < Mv%, \li,t 

6. ELi(4 + 4) <1, Vz,t 

7. v%<x\,\fi,t,k 

9. 4 = 4 = o,viGAK} 



10. 4 G {0, 1} 

IL 4e{o,i} 
12. 4e{o,i}- 

Here y% = ZU '^^Ph " ^ Zl=i TZ=i hmP. 



ON, 



q^k 



M is a large positive constant, and the constraints have 
the following interpretations: 

1) No negative power is allowed. 

2) Every node in the destination set is required to 
have decoded the data by the end of time slot T. 

3) If a node has not decoded a message by the end 
of time slot t, that node is not allowed to transmit 
that message at time t + 1. 

4) 4 is forced to be if message k is not decoded 
in time slot t. 

5) p^^ is forced to be 0, if node i is not transmitting 



message k at time t (i.e. if v-. 



0). 



6) A node cannot transmit and receive at the same 
time and can only transmit or receive a single 
message at each time slot. 

7) v^^ is forced to be 0, node i is not allowed to 
transmit message k at time t (i.e. if x| = 0). 

8) Only sources have the message at the beginning. 

9) No one else has the message at the begiiming. 
10) X, z and v are binary variables. 

We call this optimization problem MCUE, for multiflow 
cooperative unicast with Energy Accumulation. 

III. Special Case of A; = 1 

In this section we consider MCUE for the special 
case of A: = 1 and prove the problem can be solved 
optimally and in polynomial time for this special case. 
We also provide a polynomial-time algorithm to achieve 
the optimum solution. 



The problem can then be formalized as a combinato- 
rial optimization problem: 



Theorem III.l The optimal solution for MCUE is a 
simple path for k = 1, but not necessarily so for k > 2. 

Proof: The claim can be proved by induction on 
T: For delay T = 1, the claim is trivially true, as the 
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Level 1 Level 2 Level 3 Level T 

Fig. 1: Applying the multicommodity flow technique for 
unicast cast 



optimal solution is direct transmission from the source, 
s, to the given destination, d. Let us assume the claim 
is true for T = t — 1. To complete the proof, we need 
to show the claim holds for T = t. Pick any node in 
the network as the desired destination d. If the message 
can be transmitted from source s to d with minimum 
energy in a time frame less than t, then an optimal simple 
path exists by the induction assumption. So consider 
the case when it takes exactly T = t steps to turn 
on d. The system is memoryless, so d must decode by 
accumulating the energy transmitted from a set of nodes, 
V, at time t. This can be represented as ^ Pv.thdv, ^ ^■ 

We observe that there must exist a node Vo € y whose 
channel to d is equal or better than all the other nodes 
in V. Therefore, given hdv„ > hdv.yvi e v\{fo} then 
E Pvd^dv^ > Yl Pvdhdv, > ^- In other words, if we 

add the power from all nodes in v and transmit instead 
from Vo, our solution cannot be worse. Vq must have 
received the message by time t — 1, to be able to transmit 
the message to d at time t. We know by the induction 
assumption that the optimal simple path solution exists 
from source to any node to deliver the message within 
t — I time frame. Thus, for T = t, there exists a simple 
path solution between s and d, which is optimum. ■ 
Considering the above theorem, the MCUE problem 
formulation (for the special case of /c = 1) reduces to: 



min Ptotal 
s.t. 1. 
2. 

3. 

4. 
5. 
6. 
7. 



- Z^i=l l^i=l Pit 
Pit > 0, Vi, t 
XdT+l = 1 

n 

-M{1 - Xit+i) < J2 hjxPjt 

Pit < Mxit, yi,t 

Xsl = 1 

Xii = 0, Vi ^ s 

Xit G {0, 1} 



(4) 



This can be solved optimally in polynomial time using 
dynamic programming. Let C{i,t) be the minimum cost 
it takes for source node s to turn on i, possibly using 
relays, within at most t time slots. Then we can write: 



C{i,t) 



min [C{j,t 

j€Nr{i) 



1) + Wji] 



(5) 



where 



with C{s,t) = 0, for all t and C(i, 1) = Ws 
Nr{i) is the set that contains i and its neighboring nodes 
that have a non-zero channel to i, Wji represents the 
power it takes for j to turn on i using direct transmission. 
Thus the solution to (4) is given by C(d, T) and its 
computation incurs a running time of O(n^). 

IV. INAPPROXIMABILITY RESULTS 

For A; = 1, we proved in Theorem III, that the optimal 
solution is a simple path. For > 2, we can consider 
the following counter-example to argue that the solution 
is not necessarily a single -path. Consider the scenario 
shown in Figure 2, where T = 3, where the edge 
weights are equal and the edges shown in gray show 
strong interference. The red nodes cannot by themselves 
transmit the message to d2, as it causes interference for 
di and ^3 preventing them from being able to decode 
the data. However, they can cooperate with each other, 
by each sending with half power to get the message to 
d2 without causing too much interference for the other 
destinations. 




Fig. 2: An example of fc > 2, with T = 3, where the optimal 
solution is not a single path. 



To investigate the complexity of MCUE, let us start 
by looking at a sub-problem. Imagine a one hop setting 
of k source nodes and their corresponding k destination 
nodes, with no relay nodes. Due to interference, not 
all sources can transmit simultaneously. The task is 
to schedule the sources appropriately, so that everyone 
can get their message delivered to their corresponding 
destination within a time delay T. The problem is to 
find the minimum such T. Let us call this problem 
MOSP, for multi-source one-hop scheduling problem^. 

^This is essentially the problem considered in [21], though no proof 
of complexity is given in that paper. 
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It is important to note that MCUE is at least as hard as 
MOSP. Thus, any hardness results obtained for MOSP 
imply hardness of MCUE. 

In this section, we derive inapproximablity results for 
MOSP by showing that any instance of minimum graph 
coloring problem [19] can be reduced to an instance of 
MOSP 

Lemma IV.l MOSP is o(n^/^~'^) inapproximable, for 
any e > 0. 

Proof: Given an instance G{V,E), |F| = n, of the 
minimum graph coloring, we construct a bipartite graph 
G' , with the bi-partition X and Y with |X[ = |y| = n. 
For each node vi G G, we place two nodes Ui ^ X and 
u[&Y and connect them with an edge {ui,u'j). Also for 
every edge in G, Cij = {vi,Vj}, place two edges {ui,Uj) 
and {uj,u'j) in G'. We assign Ui and u[ to be a source 
and destination pair respectively for all i. We set equal 
edge weights for all the edges in G' and set > 1 to 
get an instance of MOSP 

A simple example is shown in Figure 3. Notice that 
the gray edges in the figure represent interference, and 
by setting 6* > 1, a message can be successfully decoded 
if and only if there is no interference at that node. 




Fig. 3: Example construction of G', for a given G. 



This in turn means two sources in G' can simultaneously 
transmit if and only if there is no edge in between 
them in G. Thus, the set of nodes that are transmitting 
simultaneously in G' correspond to an independent set in 
G. Consequently, the optimal solution to MOSP is equal 
to the minimum graph coloring of G, which is known 
to be o(n^/^~'^) inapproximable [20]. The following 
theorem follows by noticing that MOSP is a special case 
of MCUE. 

Theorem IV.l MCUE is o(n^/''^^) inapproximablejor 
any e > 0. 

Notice that the inapproximability result, given by 
Theorem IV.l, is stronger than, and implies, the NP- 
hardness result. In other words, it implies that not only 



finding the optimal solution is NP-hard but finding a 
polynomial time approximation algorithm that approxi- 
mates the optimal solution to MCUE with a factor of 
o(n^/^~'^) is also NP-hard. 

V. Performance Bounds 

In section IV, we proved that MCUE problem is 
in general inapproximable. However, it was shown in 
section III that the problem can be solved optimally and 
in polynomial time for the special case of A; = 1. In 
this section, we use the results of section III to obtain 
performance bounds for MCUE. 

A. An Analytical Lower Bound 

In this section we establish a lower bound on the 
optimum solution to MCUE. 

To get a better intuition for this lower bound, let us 
start off by considering the optimal solution to MCUE 
for the case when there is only one flow present in the 
network. As before, we have n nodes and a channel 
H, but this time the source s wants to transmit the 
message to a particular destination d, using the minimum 
energy within a given delay constraint T. The system 
is cooperative in that other nodes in the network, may 
be utilized as memoryless energy accumulating relays to 
help achieve the minimum energy goal. Based on section 
III, the solution can be found by calculating G{d,T) 
where C{d,T) is defined as per equation (III). 

To find a lower bound for MCUE for a general case 
of r flows, with source-destination pairs {(sfc, (ifc)}^_^, 
all with the same delay constraint T, we notice that 
the cost paid by optimal MCUE to cover each node 
cannot be lower than the optimal minimum cost paid 
by each source to cover its corresponding destination 
dk in the absence of other interfering flows. Based on 
that observation we derive the following lower-bound, 
LB(T), for the OPTmcue for r flows when the delay 
constraint is T: 

r 

LS(r) = 5]c(4,r) (6) 

k=l 

where G{dk,T) is defined as per equation (III). In 
other words, C{dk,T) calculates the minimum cost of 
optimal single flow transmission to cover a destination 
dk, starting from its corresponding source under a delay 
constraint T. LB(T) takes the sum of those costs and 
use it as lower-bound - since we know OPTmcue has 
to cover all these flows and cannot do so any better than 
the optimal solution for a single flow. Computing this 
lower bound incurs a running time of O(n^). 
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B. An Analytical Upper Bound 

In this section we establish an upper bound on the 
optimum solution to MCUE, for the general case of r 
flows, with T >r. 

The upper bound is established by considering the 
multiplexing solution. At the extreme end of T = r, 
we would allow one time slot for each of the r flow 
to transmit its message, while the other flows are silent. 
For a general time T(> r) we break the time into r 
blocks T = (ri, Tfc, ...r,.), such that Yl\=i'''k = T. We 
assign each block to one of the flows, while the other 
flows are silent. We calculate C{dk,Tk), defined as per 
equation (m). For a given tuple T, the summation of 
the total energy required by all flows to complete their 
transmission can be achieved by calculating: 

r 

UB{T) = J2C{dk,rk) (7) 
fe=i 

This sum would provide an upper bound for OPTmcue- 
For a general T < r, v/e will have (^Zl) possibilities 
for assigning the time slots to different flows. The upper 
bound is calculated as follows: 

UB{T) = mm UB{T) (8) 

To compute this upper bound we need to carry on 
the computation for calculating a single flow MCUE, 
discussed in section m, (^Zl) x f times. Thus the upper 
bound incurs a miming time of 0{nP). 

VI. A Polynomial-time Heuristic 

In this section we propose a polynomial time heuristic 
for MCUE, the performance of which is later evaluated 
against that of the bounds established in section V. We 
call this polynomial time heuristic MCUH, for multiflow 
cooperative unicast heuristic. 

To recap, consider a network, G, with a total of 
n nodes, / = {l,..,n}. Assume we have r source 
nodes, labeled S = {si, S2, Sr}, and r corre- 
sponding destination nodes, V = {di,d2, ■■■,dr}. The 
source-destination nodes can be thought of as pairs, 
{{sk, dk)}k=i, all with the same delay constraint T. The 
goal is to deliver a unicast message from each source to 
its corresponding destination, possibly using other nodes 
in the network as relays. The objective is to do so using 
the minimum amount of sum transmit power and within 
a given delay constraint. 

The MCUH algorithm works greedily by scheduling 
flows one by one. Each flow is given more slots than 
its previous flows, to ensure a feasible solution always 



exist. That means the algorithm works for T > r. 

Each flow, with the exception of the final flow, uses 
more power than required to deliver its message. This 
is achieved by assigning a higher threshold to that 
flow when scheduling the flow. After scheduling, the 
nodes that will be transmitting at each time slot and 
the power they use for transmission is passed on to the 
next flow. Each flow, when scheduling itself, will ensure 
that its transmission will not disturb the transmission of 
previously scheduled flows. A lower threshold is used to 
check for disturbance, than the one used for scheduling 
the flow itself. Let us now look at the details of the 
algorithm. 

We schedule the r flows greedily, starting from the 
one that causes the least disturbance. Without loss of 
generality, let us assume that we are scheduling the flows 
in the order 1 to r. All flows need to be scheduled within 
a total of T time slots. For flow 1, we assign a delay 
constraint of Ti time slots for transmission, for flow 2, 
we assign T2 and so forth, such that: 

1 <Ti <T2 < ... <Tr = T (9) 

Recall that time-slots are in unit durations, thus can 
increment in integer units. Therefore, each flow has at 
least one more time slot at its disposal than its immediate 
predecessor, ensuring that a feasible solution always 
exists. This also means that the algorithm works for 
T>r. 

In section II, we defined 9 to be the decoding threshold 
as per equation (1). For this multi-flow setting, each flow 
is assigned its own 9 value, such that: 

9i>92> ... >9r = 9 (10) 

Flow 1 is scheduled with Ti and 9i, as per algorithm 
in section m. We store the nodes that are scheduled 
to transmit in each time slot, and their transmit power 
and their corresponding receivers in a black list B. 
According to this definition, B{t) gives us the set of 
already scheduled nodes that are transmitting at time t 
and their corresponding powers, and their corresponding 
receiving nodes. 

For the kth flow, we use a modified version of the 
DMECT-go algorithm discussed in [8], [18]. DMECT-go 
is a polynomial time algorithm that uses a deterministic 
dynamic program to optimally solve the problem of 
joint scheduling and power allocation, for any given 
ordering, in a single-flow network. In [8], [18], it was 
shown that Dijkstra's shortest path algorithm provides a 
good heuristic for ordering in a network with uniform 
distribution of nodes. 

As mentioned, DMECT-go is for a single source 
problem and does not take interference or deliverance of 
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multiple messages from multiple sources into account. 

The algorithm proposed in this section, MCUH, is a 
modified version of DMECT-go. In MCUH, we broke 
the NP hard problem of MCUE into three subproblems 
namely, ordering, scheduling and power allocation. Or- 
dering, for a vector of n nodes, is defined as an array 
of indices from 1 to n; any node that has decoded the 
message will only be allowed to retransmit when all 
nodes with smaller index have also decoded the message 
(and are thus allowed to take part in transmission). Given 
ordering, what remains to be determined is scheduling, 
and power allocation. In other words, what remains is 
deciding which nodes should take part in transmission 
of each flow, how much power they should transmit with 
and at what time slots, such that minimum energy is con- 
sumed while delay constraints are satisfied. We also have 
to ensure the interference is taken into account when 
scheduling different flows and allocating powers. The 
MCUH algorithm solves the joint problem scheduling 
and power allocation, for flow A; > 1, as follows: 

Ordering: 

Pick a subset of nodes (as potential relays) and assign 
an ordering to those nodes. Let us call this ordered 
subset Ifc = (lfc,2jt, ...,jk, ■■■,'nk), where 1/j corresponds 
to Sfc and Uk corresponds to dk- This set could for 
instance be obtained by picking the nodes that would 
have been picked if we were to run the single-flow 
algorithm of section III for flow k. In [8], [18], it was 
shown that Dijkstra's shortest path algorithm provides a 
good heuristic for ordering in a network with uniform 
distribution of nodes. 

Notice that ordering does not dictate the time slot at 
which a node should transmit. It only states that any 
node that has decoded the message will only be allowed 
to retransmit when all nodes with smaller index have 
also decoded the message. Being allowed to transmit, 
also does not indicate that a node will in fact transmit. 

Given this ordering, we now need to solve the joint 
scheduling and power allocation problem. 

Power Allocation: 

Let us call the power allocation algorithm, PAM 
for power allocation multi-flow. PAM, to be specified 
shortly, calculates the instantaneous optimal power allo- 
cation for flow k at time t, given the set of instantaneous 
senders and receivers for flow k and the set B{t) (of 
senders and receivers of flows 1 to A; — 1 and their 
corresponding powers at that time slot). The design of 
PAM is as follows: 



PAM is defined for flow k, to take as an input a set 
of transmitters (^'fe(t) = {l...ife}) and a set of receivers 
(Tlk{t) = {ik + and the set of already scheduled 

nodes for that time-slot and their corresponding powers 
13 (t), the channel between the nodes and the receiving 
threshold 9^. The objective of PAM is to minimize 
the total sum power used to transmit flow k, by all the 
nodes transmitting that flow at time t. In other words, the 
objective is to minimize X]qG*fc(t) ^'gt' while satisfying 
a number of conditions. The output of this algorithm 
is the set of for q € ^fc(i) and also the sum of 
powers, r2^(fc) = J2q£<iik{t)Pqt' which is the objective 
of the optimization. We also define a corresponding 
function of the same name PAM, that returns as 
its output for the corresponding input. In other words 

PAM{^kit),nk{t),ek,H,B{t)) = niit). pam is 

described in Algorithm 1 . 

Algorithm 1 Power Allocation Multiflow (PAM) {k,t) 

1: INPUT: ^k{t) and TZkit) for a given flow k at time 

t, B{t), H, 9k 
2: OUTPUT: p'^ for q G \E'jt(t) and the objective value 

3: Begin: 

4: if g G B{t), \/q G 7^fc(^) then 

5: ni{t) := OO. 

6: return infeasible. 
7: else 



niit)=mm Pqt s.t. 

qe^k{t) 

v/e{i,...,fe-i} 

3. Yl hvzpit-Of I ^+ E 



ue*g(t) 



>o, 



Vae{i,...,fc} 

V/€{l,...,fe-l} 



4. P% = 0, Vg G B{t). 



10: 
11 



return J^*(i), {p%}q^^^(jt). 



end if 
End 



The PAM algorithm works by first ensuring that a 
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node cannot receive a message for flow k at time t, if 
it has already been scheduled to participate in another 
flow in that time slot. This renders the power allocation 
task infeasible with the given set of transmitting and 
receiving nodes. In other words, it returns infinity and 
states the result to be infeasible. If that is not the case, 
PAM proceeds to calculate ^%{t) by solving a linear 
optimization problem. In the optimization formulation, 
constraint 1 ensures that there are no negative powers. 
Constraint 2 ensures that the nodes assigned to receive 
flow k at time t will in fact accumulate enough energy 
to decode the message, despite the existing interference. 
Constraint 3 ensures that the power being assigned 
to nodes in flow k, is not disturbing the previously 
scheduled flows. Constraint 4 ensures that a node cannot 
transmit a message for flow k at time t, if it has already 
been scheduled to participate in another flow in that time 
slot. 

Joint Scheduling and Power Allocation: 

Given this power allocation algorithm for flow k, and 
the ordering, all that remains to be done is scheduling. 
In other words, determining the set of transmitting and 
receiving nodes at each time slot that need to be passed 
on to the power allocation algorithm. 

Given the ordered subset of = 
(lfc,2fc,...,jfc,...,nfc), we define Ck{jk,t) to be 
the minimum energy needed for flow k to deliver the 
message to all the nodes up to node in t steps or 
less. We then use the following deterministic dynamic 
program to solve the joint scheduling and power 
allocation problems, optimally: 



also look at the effect of channel degradation in the 
overall performance. 

We consider a network of 100 nodes uniformly dis- 
tributed on a 20 by 20 square surface. The channels 
between all nodes are static, with independent and ex- 
ponentially distributed channel gains (corresponding to 
Rayleigh fading), where hij denotes the channel gain 
between node i and j. The mean value of the channel 
between two nodes, hij, is chosen to decay with the 
distance between the nodes, so that hij = d^j^, with dij 
being the distance between nodes i and j and rj being 
the path loss exponent. The corresponding distribution 
for the channel gains is then given by 



1 / hij (k) 



Notice that the minimum power calculated by different 
algorithms, shown on the y-axes of the graphs in this 
section, are normalized by value of 9 (rendering it unit- 
less). 



Lowerbound 

- Heuristic 

- Upperbound 



min [Cfc(zfe,t - 1) + ... 
«fcG(i,..,ifc) 



Ck{jk,t) 

+ PAM{{l...ik}, {ik + l...jfc}, Ok,H, Bit))] (12) 



where CkiikA) = PAM{lk,{'^k.:ik},Ok, H,l3{t))] 
yik G and Cfc(lfc,t) = Vt. Thus, for flow k, 

the total minimum cost for covering nodes by time 
Tk can be found by calculating Ck{nk,Tk). 

After each flow is scheduled, we save the set of 
scheduled transmitting nodes for each time-slot t and 
their corresponding powers in B{t). The algorithm is 
repeated for each flow k > 1. 

VII. Performance Evaluation 

In this section we compare the performance of the 
proposed heuristic against the analytical bounds for an 
example network with arbitrarily chosen three flows. We 



5.5 6 6.5 7 7.5 
Delay Constraint (T) 



Fig. 4: Performance of the heuristic against the analytical 
upper and lower bound. 

Figure 4 shows the performance of the heuristic 
against that of the analytical bounds. As can be seen 
the heuristic is performing close to the lower bound. 
Notice that the lower bound is an unachievable lower 
bound, in that it assumes no interference is present. This 
means that its performance is not achievable by any 
algorithm. This is more emphasized when we have fewer 
time slots available, and thus we need to use more power 
to transmit the message creating a lot of interference that 
is ignored by the lower bound. As we get more time-slots 
available to us, the performance of the heuristic and the 
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Fig. 5: Effect of channel degradation on the total energy 
consumed. 



bounds seem to converge, which is what we expect as 
the solution goes to a multiplexing solution in all cases. 

We see the effect of poor channel conditions in 
Figure 5. As expected the performance is degraded as 
the channel conditions become poor, this highlights the 
importance of having smart algorithms to minimize the 
energy consumption in such scenarios. 

VIII. Conclusion 

In this paper we formulated the problem of minimum 
energy cooperative transmission in a delay constrained 
multiflow multihop wireless network, as a combinatorial 
optimization problem, for a general setting of fc-flows 
and formally proved that the problem is not only NP- 
hard but it is o(n^/^^'^) inapproxmiable. We proved in- 
approximability by reduction from the classic minimum 
graph coloring problem. To our knowledge, the results in 
this work provide the first such inapproxmiablity proof in 
the context of multiflow cooperative wireless networks. 

We further proved that for the special case of A; = 1, 
the solution is a simple path and devised an optimal 
polynomial time algorithm for joint routing, scheduling 
and power control. We then used this algorithm to estab- 
lish analytical upper and lower bounds for the optimal 
performance for the general case of k flows, where the 
delay constraint is at least equal to k. Furthermore, we 
proposed a polynomial time heuristic for calculating 
the solution for the general case and evaluated the 
performance of this heuristic under different channel 
conditions and against the analytical upper and lower 
bounds. 
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